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Recently a real time picture of quantum field theory has been developed which allows one 
to look into the time evolution of a scattering process. We discuss two pictures for discussing 
non-equilibrium processes, namely the Schrodinger Picture and the Heisenberg picture and 
show that a Time Dependent Variational Method is equivalent to the leading order in Large-N 
approximation in the Heisenberg picture. We then discuss the dynamics of a non-equilibrium 
chiral phase transition in mean field theory in the 0(4) sigma model. We show how the pion 
spectrum can be enhanced at low momentum because of non-equilibrium effects. We then 
show how to use Schwinger's CTP formalism to calculate the inclusive dilepton spectrum 
coming from the pion plasma. We find that a noticable enhancement occurs in this spectrum, 
but that there are large numerical uncertainties due to errors connected with the finite times 
used to do our numerical simulations. 
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I. INTRODUCTION 

In the early 70's, the study of inclusive Hadronic In- 
teractions was at the forefront of theoretical research be- 
cause of the opening of two major accelerators-FNAL 
in the United States and the ISR at CERN. As data 
poured in many of us tried to understand from various 
approaches such as Regge Poles, Fireballs, and Hydro- 
dynamics how to determine the single particle inclusive 
spectrum of particles. In that atmosphere, in 1974, Pe- 
ter Carruthers founded the Particle Physics group at 
Los Alamos which initially included Dick Slansky, Geoff 
West, David Campbell, David Sharp, Mitch Feigenbaum 
and myself, all of whom were working on this problem 
when we first arrived. At that time QCD was not yet a 
well developed subject, computers were relatively archaic 
and the idea of a first principles approach to understand- 
ing what happens after two particles collide was just a 
dream. Recently this dream has gotten closer to reality, 
and as a tribute to Dick and our early efforts I would 
like to review progress we have made in this direction. I 
will begin by first reviewing two pictures for discussing 
non-equilibrium processes, namely the Schrodinger and 
Heisenberg Pictures. In the Schrodinger Picture one can 
reduce the number of degrees of freedom by using a time- 
dependent variational method with a trial density matrix 
which is Gaussian. In the Heisenberg picture one can use 
the large N expansion of the Path Integral in the Closed 



Time Path Formalism to give us a controllable expan- 
sion about the mean-field approximation. We will then 
review our model for a non-equilibrium phase transition, 
namely the 0(4) linear sigma model, and set up the time 
evolution equations for the lowest order in Large N ap- 
proximation. In this model we will obtain natural cooling 
(quenching) of the plasma from the unbroken phase to 
the broken symmetry phase as a result of expansion into 
the vacuum with boost invariant kinematics imposed. We 
find that there is a growth of unstable modes when the 
order parameter (effective pion mass) goes negative dur- 
ing the evolution. This causes a distortion of the sin- 
gle particle distribution of pions. We also reconstruct 
typical classical configurations by sampling the quantum 
density matrix and see domains. We then obtain the 
dilepton spectrum from this time evolving plasma using 
Schwinger's closed time path formalism. 



II. STRATEGIES FOR STUDYING TIME 
EVOLUTION PROBLEMS IN \<p 4 FIELD 
THEORY 

A. Schrodinger Picture 

In quantum mechanics, time evolution problems are 
usually discussed in the Schrodinger Picture. The initial 
state in the x representation is the wave function 



*(x,i) =< a;]* > 



(2.1) 
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which evolves in time according to the Schrodinger equa- 
tion: 
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(2.2) 



< <£|*„ >= ip v [ip,t] 



(2.9) 



This equation can be obtained from Dirac's Variational 
Principle Define the action 



r = y df < V\id/&t - H\V > 



(2.3) 



Minimimizing the Action then yield the Schrodinger 
equation: 



ST = -> {i— - ff}|* >= 



(2.4) 



Here one thinks of _ff as an operator which in the x rep- 
resentation one has 



P 



dx 



(2.5) 



Gaussian initial conditions for the wave function would 
be: 

^(x, 0) = exp[— a(x — xo) 2 }- 

In field theory, the wave function is replaced by a wave 
functional. For example a Gaussian Wave functional is 
given by: 

< <p\9 >=ip[(p,t] = 

exp[[- f [<p{x) - y{x)\ [ ° 1{ ^ V) - iX(x, y)] [<p(y) - 

J x,y 



+i7t(x)[ip(x) — (f(x)}] 



(2.6) 



The time evolution is now given by the functional 
Schrodinger equation : 

H = J d 3 x[-±6 2 /5^ 2 + \{V V ) 2 + V{y)] (2.7) 
Dirac's variational principle 

r = J dt <ty\id/dt- H\V > (2.8) 



with ip v given by Eq^^) The variational parameters 
have the following meaning: 

tp(x,t) = < VvlflVv >;n(x,t) =< -i5/6<p\V v > 
G(x, y, t) =< V v \<p{x)y{y)\^ v > -<p(x, t)<p(y, t) (2.10) 

The Effective Action for the variational parameters is 
r(<£,7r,G,£) = J dt <V v \id/dt- H\f v > 
= I dtdx[n(x, t)if(x, t) + J dtdxdyT,(x,y)G(x,y,t) 
dt<H > (2.11) 
where 

<H> = J dx{n 2 /2 + 2EGE + G _1 /8 + 1/2(V<^) 2 

-1/2V 2 G + l/2V"[cp]G + l/8V""[ip]G 2 }. (2.12) 

Equations of motion that result from varying the ac- 
tion are: 

w(x,t)= V 2 ip-d <V> /dip; 

<f(x, t)= 7T 

G(x,y,t)=2 J dz[E(x,z)G(z,y) + G(x,z)E(z,y)] 

±(x,y,t)=-2 J dz[i:(x,z)Z(z,y) + G- 2 /8 

+ [\v 2 x -d<V> /dG]8\x~y) (2.13) 

In tfi 4 theory if there is translational invariance, we can 
simplify these equations by Fourier transforming them in 
three dimensional spase to obtain: 



5T = 



{i--H}\V>=0 



is the starting point for thinning the degrees of freedom. 
Instead of putting the functional Schrodinger equation on 
the computer, one assumes that the wave functional stays 
Gaussian, and uses the variational principle to determine 
the time evolution equations for the one and two point 
functions involved in the description of the Gaussian. 
This approximation is equivalent to mean field theory 
and at large-N becomes equivalent to the leading order 
large-TV approximation. That is we assume a Gaussian 
trial wave functional: 



2G(k,t)G(k,t) - G z (k,t)+'ir(k,t)G z (k,t) -1 = 
T(k, t) = k 2 + m 2 (t);m 2 (t) = -/i 2 + - A / [dk]G(k, t) 

(2.14) 

One can also linearize these equations by recognizing that 
the mean field approximation for the homogeneous prob- 
lem is equivalent to a field theory with a time dependent 
mass which is self consistently determined. That is if we 
assume a quantum field obeying: 



(□ +m 2 {t))(j)[x,t) = 



(2.15) 



then one can satisfiy the equation for G(x, y, t) by choos- 
ing: 



G{x,y;t) =< ct>(x,t)4>(y,t) > 



(2.16) 
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B. Heisenberg Picture— Path Integral Approach- 
Closed time-path approach of J. Schwinger 

In the Heisenberg Picture, the operators are time de- 
pendent, and the expectation value of the fields in an ini- 
tial state are the infinite number of c-number variables. 
The infinite heirarchy of coupled Green's functions need 
to be truncated by an approximation scheme. The large 
N approximation orders the connected Green's functions 
in powers of 1/N, with the connected four point function 
going as 1/N, 6 point function 1/N 2 etc. The formalism 
for preserving causality in initial value problems was in- 
vented by Schwinger jq] and was later cast in the form of 
path integrals j| . The Generating Functional for initial 
value problem Green's functions is 

Z[J+,J-,p}= J d^ out < i|T*(exp{- J iJ_tp_})\* out > x 

< *out|T(exp J iJ+(p+)\i > 

= Tr{pT* (exp{- J iJ_y>_})T(exp J iJ+(f+)} 

(2.17) 



G ab (t,t') = 



5 2 W 



SJ a (t)SJ b (t<) 



(2.23) 



This can be written as the product of an ordinary Path 
integral times a complex conjugate one or as a matrix 
Path integral. 

Z[J+,J-} = J n^Trp{expi [(,%+] + J J +v+ ) 
-m<P-] + J J-V-)]} 

a 

(2.18) 



= gtW[Ja] 



G 21 (t,t') = G>(f,0 = iTr{p y{t)Tp{t')} con , 
G 12 (t,t') = G<(i,0 = ±iTr{p W'Mt)} con 

G n (t,t') =iTr{pT[ l p(t)Tp(t')]} C on 

= 0(i, t')G> (t, t 1 ) + Q(t', t)G< (t, 

G 22 (t,t') = l Tr{p T*[ V {t)W)\} con 

= e(t', t)G> (t, t 1 ) + 9(t, t')G< (t, t') 

(2.24) 

We notice that G F = G n (t,t') and Gf* = G 22 (t,t'). 
We also will need the relationships: 

G ret (t,t')=ie(t-t')mt),$(t')]± 

= e(t-t')[G>(t,t') - G<(t,t% (2.25) 

G adv (t,t')= -iB(t' -t)[$(i),$(i')]± 

= Q(t' - t)[G<(t,t') - G>(t,t')}. (2.26) 

as well as the relations between the Green's functions: 

G ret (t, t')= G F (t, t') - G<(t, t') = -Gf* (t, t') + G> (t, if) 
G adv (t, t')= G F (t, - G> (t, = -Gf. (t, t') + G< (t, t') 

(2.27) 



1. Large-N Approximation 

The method for reducing the number of degrees of 
Freedom in the Heisenberg picture is the large- N approx- 
imation H If we have an ./V-component scalar field with 
Lagrangian: 



where < ip%, i\p\ip2, i > is the density matrix defining the 
initial state. We use the matrix notation: 

V a = ; « = 1,2 (2.19) 

with a corresponding two component source vector, 

J a = ( J T + ] ; a = 1,2. (2.20) 



A 



2 V "' ,I ' A " 87V 1 ' ' ' 



(2.28) 



On this matrix space there is an indefinite metric 

c ab = diag (+1, -1) = c ab (2.21) 

so that, for example 

J a c ab (p b = J+Lp+ - J-<p- . (2.22) 

From the Path Integral we get the following Matrix 
Green's Function: 



we can rewrite this as: 

Lci[$, x}= -7^( D + + yX (| + M 2 



where i — 1, .... N and 



X = ~M H , 

2iV 



(2.29) 



(2.30) 



If /z 2 > 0, spontaneous symmetry breaking at the clas- 
sical level. At this minimum the O(N) symmetry is spon- 
taneously broken, x — an( i there are N — 1 massless 
modes. Small oscillations in the remaining i = N (ra- 
dial) direction describe a massive mode with bare mass 
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equal to \[2\x = \f\vn. The Generating functional for all 
Graphs is given by [|6|: 



Z[j,K]= J d<t>dxexp{iS[cl>,x]+i j\j<t> + K X ]} (2.31) 
Perform the Gaussian integral over the field 4> 

Z\j,K]= f d X eMiNS eff \x,j,K]} (2.32) 



approximation boost invariant kinematics which result 
from a hydrodynamic picture where the original plasma 
is highly Lorentz contracted. In the linear sigma model 
treated in leading order in the 1 /N expansion the theory 
has a chiral phase transition at around 160 Mev and we 
choose the parameters of this theory to give a reasonable 
fit to the correct low energy scattering data. We obtain 
natural quenching for certain initial conditions as a result 
of the expansion process. 



where 

S eff = [ ^{^G- 1 [xb-+^X+^x(|+M 2 )+^TrlnG- 



A. review of the linear a model 



G- 1 \ X ](x,y) = {a+ x }6(x-y), 



(2.33) 



Because of the N in the exponent one is allowed to per- 
form the integral over x by stationary phase. This leads 
to an expansion of Z in powers of 1/N the lowest term 
(stationary phase point) is related to the previous Gaus- 
sian (Hartree) approximation. The effective action of the 
leading order is: 

S e //[<M = S cl [<t> >X ] + y TrlnG- 1 !*] . (2.34) 
Varying the action leads to the mean field equations: 

{O + x }4> = (2.35) 



IN i 

We notice that this is the same equation found in the 
Gaussian approximation with m (t) being identified with 
X- 



The Lagrangian for the 0(4) a model is: 

L = -d$ ■ <9<i> — -A($ • $ — v 2 ) 2 + Ha. (3.1) 
2 4 

The mesons form an 0(4) vector 
$ = (a,Wi) 

As we discussed earlier in our discussion of the large- N 
approximation we introduce: 

X = A($ • $ - v 2 ) 

and use the equivalent Lagrangian: 

L2 = -\u^ + x)^ + ^ x +\xv 2 +Ha (3.2) 

The leading order in 1/N effective action which we 
obtain by integrating out the </> field and keeping the 
stationary phase contribution to the x integration is 

x\=J d 4 x[L 2 (<P, X; H) + % -Ntx InGo- 1 ] (3.3) 



III. DYNAMICAL EVOLUTION OF A 
NON-EQUILIBRIUM CHIRAL PHASE 
TRANSITION 

One important question for RHIC Experiments is- can 
one produce disoriented chiral condensates (DCC's) in a 
relativistic heavy ion collision? Recently, Bjorken, Ra- 
jagopal and Wilczek and others proposed that a nonequi- 
librium chiral phase transition such as a quench might 
lead to regions of DCC's |Q. The model Rajagopal and 
Wilczek considered was the 0(4) linear sigma model in a 
tree-level approximation, where a quench was assumed. 
Two deficiencies of that model were its classical nature 
(it could not describe ir-ir scattering V, and the quench 
was put in by hand. Our approach |8|] instead was to 
look at the quantum theory in an approximation that 
captures the phase structure as well as the low energy 
pion dynamics. We also used the natural expansion of an 
expanding plasma to cool the plasma and built into our 



G^(x iy )=i[a + X (x)} S\x-y) 

This results in the equations of motion: 

[□ + X (x)]m =0 [□ + X (x)W = H, (3.4) 

and the constraint or gap equation: 

X = -Xv 2 + \(o 2 + it ■ tt) + XNGo(x, x). (3.5) 

We will introduce fluid proper time and rapidity vari- 
ables to implement the kinematic constraint of boost in- 
variance. 

To implement boost invariance we assume that mean (ex- 
pectation) values of the fields $ and x are functions of r 
only. 
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X (r) = X(-v 2 + $ 2 (r) + AT Go ( Xj x - Tj T )) , (3.6) 

To calculate the Green's function Gq(x, y; r, r') we first 
determine the auxiliary quantum field <^(x, r) 

(t- 1 ^ rd T - r~ 2 ^ - 9i + x{x)) <t>{x, t) = 0. (3.7) 



Go(x,y,T,T') =< T{<p(x,r) <t>{y,T)} > . 
We expand the quantum field in an orthonormal basis: 

<f>(r),x±,T) = J[d 3 k](exp(ikx)f k (T) a k + h.c.) 

where kx = k v T] + k±x±, [d 3 k] = dk v d 2 k^/(2ir) 3 . The 
mode functions and x obey: 

/ k + (^ + fci+x(r) + i ^)/k = 0. (3.8) 

x(t) = x(-v 2 + <5> 2 {t) + ±-N J[d\]\Mr)\ 2 (l + 2n k )). 

(3.9) 

when x goes negative, the low momentum modes with 

+ kl < \x\ 

grow exponentially These growing modes then feed back 
into the x equation and get damped. Low momementum 
growing modes lead to the possiblity of DCC's as well 
as a modification of the low momentum distribution of 
particles. To fix the parameters of this mode we use the 
PCAC relation: 

d^{x) = f n ml-K i {x)=HK i (x). (3.10) 

and the definition of the broken symmetry vacuum. 

Xo^o = Jt4<t = H 

°o = U = 92.5MeV 
ml = -Xv 2 + A/ 2 + XN [ [d 3 k] 



The mass renormalized gap equation is 
x(r)-m 2 =-A/ 2 + A$ 2 (r) 
+-N Ad 3 k]{|/ k (r)| 2 (1 + 2 n k ) 



2^ 



X is chosen to fit low energy scattering data 

We choose our initial data (at r = l)so that the system 

is in local thermal equilibrium in a comoving frame 



1 



n k 



(3.12) 



where O = 1/T and E° = 



k\ 



x(to)- 



The initial value of x is determined by the equilibrium 
gap equation for an initial temperature of 200MeV and 
is .7 fmT 2 and the initial value of a is just The phase 
transition in this model occurs at a critical temperature 
of 160MeV. 

To get into the unstable domain, we then introduce fluc- 
tuations in the time derivative of the classical field. 
For tq = lfm there is a narrow range of initial values 
that lead to the growth of instabilities .25 < \&\ < 1.3. 

The results of numerical simulations described in Q 
for the order parameter x ar e shown in figure 1. 





I lcta/dt=-l 


do/dl=0 




I 




- V 





11.0 

X (fm) 



FIG. 1. Proper time evolution of the x field for two differ- 
ent initial values of a. 

Fig.l displays the re sults o f th e numerical simulation 
for the evolution of x (|3.8| )-(^S|). We display the aux- 
iliary field x m units of /to~ 2 , the classical fields i> 
in units of fm -1 and the proper time in units of fm 
(\frn~ 1 = 197MeV) for two simulations, one with an 



instability (<r| 



T) and one without (<j\ 



0). 



(3.11) 



We notice that for both initial conditions, the system 
eventually settles down to the broken symmetry vacuum 
result as a result of the expansion. We also considered 
a radial expansion and obtained similar results Q. In 
the radial case, the outstate was reached earlier, but the 
number of oscillations where x became negative was sim- 
ilar. To determine the single particle inclusive pion spec- 
trum we go to and adiabatic basis and introduce an in- 
terpolating number operator which interpolates from the 
initial number operator to the out number operator. In- 
troduce mode functions /P which are first order in an 
adiabatic expansion of the mode equation. 
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fO 
Jk 



-iVk(r) 



rEJv 



dyk/dt = uj k , 



(3.13) 



-1/2 



]d 3 k](exp(ikx)f^(T) a k (r) + h.c. 



(3.14) 

In terms of the initial distribution of particles no(k) 
and [3 we have: 

n k (r) = f(k Vl k±,T) =< a\(T)a k (r) > 

= n (k) + |,3(fc,T)| 2 (l + 2no(*)). (3.15) 

where /3(fc,r) = - ^-fk), n k (r) is the interpo- 

lating number density. The distribution of particles is 

/(^,r) = ^ idfci ^ v (3.16) 

Changing variables from (t),^) to (z,y) at a fixed r 
we have 

irdz dx\ J f(k„,kx,T) 



E 



d 3 k ndy dk\ 

= A± dk v f(k v , fc_L,r) 



f{k rh k±,r)k fJ -da tJl . 

(3.17) 

To compare our field theory calculation with some 
standard phenomenological approach, we considered a 
hydrodynamic calculation with boost invariant kine- 
matics and determined the spectrum assuming that at 
hadronization the pions where at the breakup tempera- 
ture T — (as well as T = lAm n ), with the distribu- 
tion given by the Cooper-Frye-Schonberg Formula |[l| 

d 3 N d 3 N f 

E i«=^ = M^ < 3 - 18 > 

Here g(x, k) is the single particle relativistic phase space 
distribution function. When there is local thermal equi- 
librium of pions at a comoving temperature T c (t) one 
has 



g(x, k) = ^{exp^^/Tj - l}" 1 . 
The comparison is shown in Figures 2 and 3. 

do7di=-l initial conditions 



(3.19) 



I 1 

1 j mean field approximation 

1 i 






. T c =1.4m 





FIG. 2. Single particle transverse momentum distribu- 
tion for a = — 1 initial conditions compared to a local 
equilibrium Hydrodynamical calculation with boost in- 
variance. 

do/dx-0 initial conditions 



x mean field approximation 




FIG. 3. Single particle transverse momentum distribu- 
tion for a — initial conditions compared to a local 
equilibrium Hydrodynamical calculation with boost in- 
variance. 

We therefore find that a non-equilibrium phase transi- 
tion taking place during a time evolving quark-gluon or 
hadronic plasma can lead to an enhancement of the low 
momentum distribution of pions. 



B. Determination of the Effective Equation of State 

Equation of state is obtained in the frame where the 
energy momentum tensor is diagonal- we are already in 
that boost invariant frame : 

7j t „ = diag{e,p n ,p±} 

When we have massless goldstone pions in the a model 
(H = 0) then \ goes to zero at large times. In the spa- 
tially homogenous case: 

< Too >= e < Tij >= pSij 

The equation of state becomes p = e/3 at late times even 
though the final particle spectrum is far from thermal 
equilibrium. 



2.0 4.0 

Kn If™') 
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FIG. 4. Equation of state | as a function of r for the 
massless a model where we start from a quench. 



C. Dephasing and looking for DCC's 

As we have shown in jnj], dephasing justifies the re- 
placement of the exact Gaussian p by its diagonal ele- 
ments. At large-N or in mean fieldtheory the density 
matrix is a product of Gaussians in space: 



{Vk\PeffWk) 

(27r£)-*exp 



\2 



<Pkf 



(3.20) 



After a short while because of dephasing, he Gaussian 
distribution off the diagonal <p' k = <fk is strongly sup- 
pressed 



2fcn(fc) 



This is shown in Fig. 5. We find no support for 
"Schrodinger cat" states in which quantum interference 
effects between the two classically allowed macroscopic 
states at v and —v can be observed. 
An ensemble may be regarded as a classical probability 
distribution over classically distinct outcomes The parti- 
cle creation effects in the time dependent mean field give 



rise to strong suppression of quantum interference effects 
and mediate the quantum to classical transition of the 
ensemble. 




FIG. 5. The Gaussian p eff for k = .4 from Ref. (k]] 
illustrating the strong suppression of off-diagonal compo- 
nents due to dephasing. 

If we project the density matrix onto an adiabatic num- 
ber basis, we can reconstruct classical field configurations 
from the diagonal density by replacing the field operator 
a(k) by 



a(k) -> [(n(k)] 



l/2 e #(fe) 



with n(k) obtained by throwing dice on the density ma- 
trix and <p being randomly chosen between < 4> < 2ir. 
Typical field configurations as a function of r (averaging 
over angles) are shown in Fig. 6. 



2.0 



0.0 ^ 



-0.0 




0.0 



10.0 
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30.0 
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FIG. 6. Four typical field configurations drawn from 
the same classical distribution of probabilities. 



IV. INCLUSIVE DILEPTON PRODUCTION AND 
SCHWINGERS CLOSED TIME PATH 
FORMALISM 

Schwinger's CTP Formalism is designed to allow one 
to calculate expectation values of operators in the initial 
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density matrix. One quantity we are interested in for 
obtaining an effective hydrodynamics is the expectation 
value of the energy momentum tensor: 

< in\T^(x)\in >= (e + p)u»u v - pg» v (4.1) 

where T fJ,v (x) is the Field Theory energy momentum ten- 
sor. Also by fourier transforming this energy momentum 
tensor and looking in a comoving frame, we can ask how 
much energy is in the "free" part of various components 
and define an equivalent number of quanta by dividing 
by fius k for each specie. 

If we consider the inclusive production of electron 
positron pairs the probablility amplitude is 

< e-(k,s)e + (k',s')X\i >=< X\b™d!ffi\i > . 
The inclusive distribution function for dileptons : 
E k EL d 6 N 



Using the relations between b, d to \& and the Free "out" 
fields we obtain 

in < A J d 3 xid 3 x 2 d 3 x$d 3 X4 
e^-^{u+ s ^ out \x 2 )}{^ out)+ {x 4 )u Ks } 

e ifc' (xi -x 3 ) | w + ^(out) ( X3 )}{#(°«*+) ( Xx ) Vk , s ,}\i> 

Now using the weak asymptotic condition [H that 

^It^oo = Z 1,2 ¥ out) (4.2) 

inside of matrix elements as well as the equation of mo- 
tion of the spinors and the identity: 



f t2 dF 
J _ = F (t3 )- f((l ), 



(4.3) 



we obtain: 

ou t < e~{k,s) e + (k',s') X\P 1 P 2 > in = 
i 2 Z- x J ' d 4 Xl d 4 x 2 e l(kX2+k ' x ^ 

u k>s V out < X\T{^(x 2 )^(x 1 )}\P 1 P 2 > m Vv k > >s , (4.4) 

Squaring this amplitude and summing over X we obtain: 

E k E' k d e N 

m m [d 3 k][d 3 k'} [ ' ' ' ' 1 



in < PiP 2 \T*{y{x 3 )*{x i )}T{*{x 2 )*{x 1 )}\P 1 P2 > in 

(4.6) 

is precisely the type of Green's function that is obtained 
from the generating functional of Schwinger's CTP for- 
malism 

The Lagrangian we will use to determine this 4 point 
function is the 0(4) linear a model + Electrodynamics. 

This Lagrangian has 3 pieces: The mesons form an 
0(4) vector $ = (%i,a). This strongly interacting La- 
grangian is given by 

1 Y 2 1 

Lstrong = + X)<t>i + ^ + + H<7 (4.7) 



X = A($ • $ - v 2 ). 
To this we add the free lepton and Photon Lagrangian: 

L = --F^F^ — -—(#• A) 2 + *[i 7li d» - m]* (4.8) 
4 la 

The interaction of the Photons with the pion plasma 
and the leptons is given by 

-eV'fVAn + CA. (4.9) 

If we treat the electromagnetic interactions perturba- 
tively in e 2 and the pions in the mean field approx- 
imation we obtain the graph shown in the Figure 7 

X 3 v X 1 




'4 



■2 



2 



in < P 1 P 2 \T*{^(X 3 )^{X4)}T{^{X2)^{X 1 )}\P 1 P 2 > in X 

(4 

The matrix element involved here, 



FIG. 7. Leading contribution from the plasma to dilep- 
ton production. The 4 fermion Graph is to be evluated 
using the Matrix CTP Green's functions. 

x The inverse propagators in the LSZ representation lop 
off the external legs and put the leptons on mass shell. 
One is left with: 

■5) ^^ = ^(*^ W ^*^ < 4 - 10 > 
where 



8 



M^(k, s; feV) = v{k', s')^u{k, s)u(k, s)Yv{k', s') 

w^{k,k') = wl„ + wl„ + wl v 

= ie 4 J d%[ d*y»J d 4 Zl J dSs e «C*+*0(*-i/0 

x[D™J(y u y 2 )ir A21 (y 2 , z 2 )^f(z 2 , Zl ) 
+£>^ (yi)2/2 )n^(y 2! z 2 )^f(z 2 , Zi) 

+£>^(y 1 ,lft)n^( W ,%)i^l(22,«i)] 

If we want the invariant Mass distribution function 
when 



M 



q = k + k 



we obtain 



d,N a°d,N 
A 24^=^(9)^(g). (4.11) 



d 4 q dM 2 d 3 q 



where 



-Rm^(9)= 



(2tt 



[ ^ [ ^S\ q -k-k')L, v (k,k') 

' 2 -d-^) 1/2 (i + ^)(^-^V) 



(4.12) 



If we were doing an ordinary perturbation theory cal- 
culation analytically we could use the translational in- 
variance of the polarization tensor: 

n£„(«2,* a )= f[d 4 q}e^y^u>M (4.i3) 



and the representation of the free Photon propagator in 
Feynman Gauge: 

D Pla ,{z2,Zi)= J[d 4 k]e-^- z ^-^- (4.14) 



to obtain: 

W^(k,k')= -ie 



= -ie 4 VT J ^-5 4 (q-k-k>)n^(q 2 ) (4.15) 

The other terms have the photon on mass shell so they 
give no contribution to the particle production rates. 

In this homogeneous case we therefore obtain the usual 
result p) pi: 



q°dN 



dM 2 d 3 qVT (2tt) 6 3 



5^ ( i-^a + ^), 



(4.16) 



Let us first look at the case of a thermal plasma where 
we can calulate everything analytically. The vacuum po- 
larization graph can be found using the following expan- 
sion of the pion field to calculate the Finite Temperature 
pion Green's functions: 



{x,t) = 



[d 3 k] 
2u k 



exp(ikx)a k + exp(— ikx)bl (4.17) 



The creation and annihilation operators obey the com- 
mutation relations: 



[a k ,al] = [b k ,bU = (2ir) 3 5 3 (k~k') 



(4.18) 



And the phase space number densities and n k are 
defined by 



{a\a k ,)= (2Tr) 3 n+S 3 (k~k') 
(b{b k ,)=(2n) 3 n-5 3 (k-k') 



(4.19) 



so that the total number of positively charged particles 
is given by: 



V = / d 3 k{a\a k ) = I d 3 xd 3 kn+ (4.20) 



For the case of neutral plasma in thermal equilibrium at 
inverse temperature (3 we have 

< = n k = l/(e?» 1) 

Using the fact that for a free pion gas: 



{<tJ{x)<t>{y))t h 
d 3 k 



2w k (2TT) 3l! ' 
{<l>{x)^{y)) th 
d 3 k 



+ e ik-{x-y) + ^ + ^-^-ifc.^-y)] 



3 L % e ik < x - y ) + (1 + n+)e- ik < x - y )] , (4.21) 



2uJ k (2ir) 



one finds that the vacuum polarization tensor is given 



by 



d 3 /ci 



d 3 k 2 



r(2^) 4 



l(2tt) 3 2w 2 (2tt) 3 
{(ki - fc 2 ) M (fci - k 2 ) v 

i n kl n k 2 S4 (l~ fc l ~ fc 2) 

+ (l + ^ i )(l + n+) ( 5 4 ( (Z + fc 1 + fc 2 )] 
(fci + fc 2 ) P (fci + fc 2 )^ 

K (1 + % 2 ) + < (1 + <)]<5 4 (<Z - fci + fc 2 )} , 
: (4.22) 

The three terms in the above relation correspond to pion 
pair annihilation, creation and bremsstrahlung, respec- 
tively. The delta functions show that only the annihila- 
tion process survives for M 2 above the dilepton thresh- 
old. We therefore obtain: 







-^>(g) = <z 2 (l - ^) f + ^n+n"_ fe , (4.23) 



with w d 



(go ± q^Jl - 4m 2 /M 2 )/2 and g 



a/M 2 + g 2 . Thus the the dilepton production rate from 
a homogenous thermal plasma of pions is given by: 



1 dN, 



(2) 



a 2 B M 



VT dM<Pq 48tt 4 q 



e+e- 



(1 



M 2 j l- 



(4.24) 



-G^(k,p;t x ,t y )= A{k)K+ v (p) + K- v (k)A(p) 

-N +v {k)M+»(p) - M-^{k)N- v {p), 

(4.29) 



and 



B{k) n = Q is = 

±ikjD(k) [i = v = j 

-(±)ikiC(k) n = i v = 

kikjA(k) fi = i v = j 



(4.30) 



A. Finite Time Effects 

For our numerical simulation of the interacting plasma 
we can only follow the time evolution of the plasma for 
a fixed time 2T which is typically about 100 fermis. 

Therefore we need to make sure we have the causal 
formulation so that dilcptons at time T only get con- 
tributions from processes occuring at times t < T. So 
we should not use Feynman propagators but instead use 
the CTP matrix propagators when we investigate the 
theory at finite times. Also the interacting plasma is 
not time translationally invariant so we must use a non- 
covariant formalism. In the mean field approximation (as 
well as for the pion gas) there is factorization of four point 
functions so that the current-current correlation function 
takes the form 

< J»(x)J^(y) > 

=< &(x)<f>(y) >< d^(x)d v &{y) > 

- < &(x)d"$(y) >< d»$(x)&(y) > 

- < (x)$(y) >< <f>(x)d v &(y) > 

+ < d**&(x)d l '$(y) >< $(x)&(y) > . (4.25) 

If we insert the mode expansion of the charged pion 
fields: 



$(x,t) = J [d 3 k] \cKp(ikx)f k (t)a k + cxp(-ikx)f k r (t)b t k 

(4.26) 

and define phase space number densities N k and N k by 
(a{a k ,)= (2t:) 3 N+S 3 (k-k>) 



(6t6fc/)= (2n) :i N-5 3 (k-k') 



(4.27) 



we obtain: 

< J^{x)J v {y) >= 

i d k dp i( fc_^. ( g_ ff) 1 



J (2tt) 3 (2tt) 3 .U^K,p,l x ,l y ), 



(4.28) 



M±(fc) = 



Nt{k) 



D{k) n = 

-(±)ikiA(k) fj, = i 

C(jfe) v = 

-(±)ikjA(k) v = j 



(4.31) 



(4.32) 



and 



A(k;t x ,t y ) = (1 + N k )[r k {t x )fk{t y )]* + N k r k {t x )fk{t y ) 
B(k;t x ,t y ) = (1 + N k )[ti(t x )f k (ty)]* +N k f* k (t x )f k (t v ) 
C(k;t x ,t y ) = (1 + N k )[t k {t x )f k {t y )Y +N k .r k (t x )f k (t y ) 

D(k;t x ,t y ) = (1 + N k )[t k {t x )f k {t y )Y +N k f* k {t x )f k {t v ). 

(4.33) 

Contracting G^ v {k,p) with L^ v = q^q v - g^q 2 we 
obtain: 

L^G^(kp;t x ,t y ) = 

l(q- (p + k)) 2 + (q 2 - q- q)(p + k) ■ (p + k)]A{p)A{k) 
-iq (p + k) ■ q[A(k)D(p) - D(k)A(p) 
-A(k)C(p) + C(k)A(p)} 

+q- q\A(k)B{p) + B(k)A{p) - C(k)D{p) - D(k)C(p)} 

(4.34) 

where 

q = k — p 

At the special case where q = we obtain: 

L^G^(kk; t x ,t y ) = Aq 2 k ■ kA(k)A(k) (4.35) 



Here k • k — > — 



-4m z 



in the infinite time limit. 



B. pion gas 



where 



To see what the effects of finite T might be, let us 
look at the case where everything is analytically known, 
namely the pion gas we discussed before. So we use the 
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known values of A, B, C, D appropriate to the pion gas 
where 



fk(t) = 



First let us look at the effect of just putting a finite 
time cutoff into the McLerran formula for W. Using the 
covariant form of the photon propagator and just cutting 
off the internal integrations to run from — T to T one 
obtains: 

W™ to "(k,k')=i^ / d 3 y 2 J d 3 z 2 J T ^dy 2o j\z 2o x 
e- i ^-^U^(y 2 ,z 2 ) 

4 pT pT 

= iV— I dy 2o I dz 2o e- iqo ^o-^o) x 



(4.36) 



<1 J-T J-T 

[dK}G^(K,K-q;y 2o ,z 2o ) 



Again looking only at the place q = and keeping only 
the pion annihilation contribution: 

L^G^(kk;y 2o ,z 2o ) = % 2 k ■ fc e *«*(«*o-*o) (4.37) 



we can perform perform the integration over T and ob- 
tain the factor: 



4( 



sinfco - 2uo k ]T 2 
go - 2uj k 



which in the limit T — > 00 becomes the factor 

2T x 2it5(q Q - 2uj k ) 

Thus we see that in the infinite time limit, the back to 
back dileptons get contributions only from pion pairs in 
the plasma with zero combined three momentum, each 
carrying energy — qo/2. In what follows we want to see 
how doing a causal calculation changes the way in which 
we go on mass-shell from this simple replacement of the 
delta function by a representation of the delta function. 
We will find extra terms which only slowly go to zero 
with the time T. 

We now use the 3 dimensional form for the propaga- 
tors: 

D*?(x,y) = iO(x Q - y Q ) 9tlv e-< x °-y^ x 

/ [dq] jg.(x-y) \ -iqa(xg-ya) _ p iqo(x -yo)] 
2qo [C 1 

= iQ(x - y a )e^(x«-y«) g ^ J [dq}e^< s -^A[q (x a - y )} 

(4.38) 

In our numerical simulations we assumed spatial ho- 
mogeneity so that the vacuum polarization has the form: 



U>^(xy) = i(J^(x)J v (y)) 

= J[dK] JldP^K-^-^G^lK^P^xo^o] 

(4.39) 

Inserting these into the expression for W^ v {k, k') we then 
obtain 

/T pT pyio pzia 

dy w / dz w / dy 2a / dz 2a 
-T J-T J —T J —T 

J ^_^-j e -e(2/io-2/2o) e -e(«io-«2o) e iqo(«io-2/io) 

A[\q\(y 10 - y 20 )}G^(K, K - q:y 20 , z 20 )A[\q\(z 20 - z w )} 

(4.40) 



where 



q = k + k' 



It is this contribution which persists when the time cutoff 
T —> OQ. 

For the other two contributions we obtain: 



W^ 2+3) (k,k') = -2e 4 V 

/T pT pzw pZ2o p 

dy w / dz w / dz 2a / dy 20 / [dK\l 
-T J-T J —T J-T J 



me 



iqo(zio-yio) 



[G» V (K, K - q- y 20 , z 20 ) - G» v {-K, q- K; z 2Q ,y 2Q ) 



-i\q\{vw-V2o) 



-A[|gl(z 20 -yio)]] 



(4.41) 



Things simplify dramatically at the place where q = 0. 
At that point the second and third contributions vanish 
and we have 



1 dN Aql 2Bn 



VT d n+1 q (2tt) 2 " 3 



0„ J k n+1 dk F[k,q = 



9oJ 

g4 pi pi p'yia rzio 

F[k,q = 0,q o ]=i— / dy w I dz w I dy 2Q I dz 20 

1 J-T J-T J-T J-T 



e iqo (zio - 2 ' lo) Ao(yio-y2o)[^(fc,y2o^2o)] 2 Ao(z 20 -z 1 o) 



(4.42) 



where 



A(k,t,t') = (l + n k ) 



e -iu> k (t-t') e iw k (t-t') 

+ n k - 



2u> k 



2u> k 



and 



. , .smm~(t — t') 

A {t-t') = -i- 



r«/2 



■ n — 2 



r(n/2) 



where n is the number of spatial dimensions. 



In the limit m~, — > we have 



A (t - t 1 ) = -i(t - t')e 



_ — el* — *' I 
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The result for the annihilation part for massless pho- 
tons is: 

e 4 N 2 (k 2 ) 
Knn ~ 64^( go -2^)2T X 

[16u>i + qt- 8tu 2 q 2 cos[2(2w fe - q )T] (4.43) 



with 



2 , 2 , 2 

ui k = k + to 



We can rewrite this as: 
F = — x 

ann A„2, ,2 A 



{ 



sin 2 [(2^- g o)r] ( go + 2^ fc ) 2 
( go -2o; fe ) 2 T + IQqfclT 



} 



(4.44) 



We would like to compare the representation of the 
delta function squared found here (renormalizcd to one 
at the delta function) to the result of just naively putting 
a cutoff into the covariant calculation which gave: 



sin 2 ({ (Z o-2w fc }r) 



(4.45) 



{qo - 2^ fe }2 T 2 

The CTP formalism which preserves causality instead 
gives: 

sin 2 [(2u fc - g0 )r] (q Q + 2cu k ) 2 
{ ( 90 - 2c fe ) 2 T 2 + 16 g 2 c 2 T2 1 

The last term makes a very small difference at large T. 
In obtaining this result we assumed that 

sT — > oo as e — > 0. 

For the brehmstrahlung contribution, we find 



rehm 



= -4 



e 4 n k (l + n k ) 



ttT (m 7 2 + e 2 ) 2 w 2 (q 2 - m 2 ) 2 



(-2m 7 4 + 2m 1 2 5o 2 
-2e 4 + 2 (to 7 2 +e 2 ) 



<7o 4 - 4 to 2 e 2 - 2 g 2 s 2 



TO. 



g 2 +£ 2 ) cos(2 7() T) 



+4g £ (to 7 2 +e 2 ) sin(2q T)) 



(4.46) 



so that the renormalized delta function squared for this 
case becomes: 

1 



f 

(-2 m 7 4 + 2 to 7 2 q z - g 4 - 4 to./ - 2 g 2 £ ; 



4q 2 T 2 (to 7 2 + e 2 ) 

4 1 o™ 2 „ 2 



-2e 4 + 2 (r 



2 e 2 

y fc 

2 1 „2 



2 _2 



q 2 +e 2 ) cos(2 goT) 



+4g £ (to 7 2 +£ 2 ) sin(2<?oT)) 



(4.47) 



Here we have kept both a small photon mass (to regu- 
late the infrared) as well as e. For the brehmstrahlung 
contribution one can not set e to zero. 



At m, 



this becomes: 



sin 2 (g T) 
1 2 T 2 



[1 



1 2 o 



1 



4e 4 T 2 £ 2 T 2 



(4.48) 



In order to obtain our previous results we used eT— > 
00. We now see to also have the unwanted terms going 
to zero we also need e 4 T— > 00. The actual rate of 
production of dileptons gets from this expression a con- 
tribution which goes to zero as a constant divided by the 
total time. This constant is about the size of the entire 
dilcpton production rate at a T = 20 fermis. In the pion 
gas case when T = 100 fermis, we need (in inverse units) 
e = 1 for the brehmstrahlung contribution to be not that 
big a contamination, and for e to be small enough for the 
annhilation cross section to be reasonably accurate. To 
take the limit E^Owe choose for T > T = 100: 



e- > (T/T ) 



< 6 « 1, 



in order to smoothly go to the covariant cutoff result that 
is 



sin 2 ( g0 T) 
V 2 T 2 



(4.49) 



With this form for the s, the finite qo dependent contri- 
bution to the cross section goes to zero as . Thus 
in doing numerical simulations, we find that in order to 
avoid contamination from brehmstrahlung processes we 
need to go to quite large hadronic time scales T > 1000 to 
be in the asymptotic regime for the production of dilep- 
tons which is an electromagnetic process. 

For the brehmstrahlung process, the effective S func- 
tion is independent of k so that one can do the integration 
over k for any qo to obtain: 



dk +1 (l + nfc)nfc 

K n 



(4.50) 



(2tt)' „ k 

For /T 1 = m x we find C\ = .127718 , and C 3 = .0792387 

For the creation contribution we get a result similar to 
the annihilation but with qo — > — qo, that is: 

e ' N k 



creation 



{ 



sin 2 [(2a> fe + q a )T] (-q + 2uj k ) 



+ 



-} (4.51) 



(qo + 2co fe ) 2 T 2 ' 16^ 2 T 2 
so that the renormalized square of the delta function is: 
.sm 2 [(2aJ k + q )T] (q Q -2LU k ) 2 



+ 



(qo + 2u; k yT 16«ga#T 
vs. the cutoff McLerran formula result: 

sin 2 ({q + 2uj k }T) 
(q + 2 Wfe ) 2 T 2 



-} 



(4.52) 
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Using these formula, we have evaluated the dilepton 
production at fixed time T for both a pion gas and for 
the interacting field thory described by the a model de- 
scribed above. When the effective pion mass goes nega- 
tive, there is significant enhancement of the signal, how- 
ever one can see the finite time effects are still not con- 
trolled in our present simulations. In figure 8 we show 
the finite time effects for a pure pion gas where one can 
determine the infinited time limit analytically. The time 
here is 100 fermis. In figure 9 we see that quench condi- 
tions enhances significantly the production of low mass 
dileptons over what one would find for a pure pion gas. 



10" 
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> 

"D 



i 10" 11 
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Numerical, 
Numerical, 
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0.7 0.8 



FIG. 9. Comparison of quench (negative effective mass 

regime) vs. non-quench initial conditions for an inter- 
acting field theory described by the sigma model. The 
large rise for small effective dilepton mass is the artifact 
of the finite time for the simulation. This effect slowly 
goes away as 



V. CONCLUSIONS 



FIG. 8. Finite Time effects for a pure pion gas. Here 
t f = 100/ 



We have shown how to use the CTP formalism to cal- 
culate the dilepton spectrum arising from a time evolving 
or a thermal plasma. For a plasma undergoing a chiral 
phase transition we expect a strong signal for existence 
of DCC-states in e + e~ - channel 



M inv ~ 2m n ; q±_ < 300 MeV 



This would be visible by CERES if 

kf* = 60 MeV 

In our calculations we have ignored possible important 
effects of direct two body scattering in the plasma which 
arise only in next order in the 1/N approach. A similar 
enhancement seen by D. Boyanovsky et. al. [|l5| in pho- 
ton spectrum. Another problem for us is that our result 
is influenced by large finite Time corrections which are 
apparent for the free pion gas. 
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